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1 Introduction 



The Tomonaga-Luttinger model ^ElEl of one-dimensional (ID) strongly corre- 
lated electrons gives a striking example of non- Fermi- liquid behaviour 01 |S1 HI ■ 
Since the seminal paper by Haldane (6j, where the notion of the Luttinger liq- 
uid (LL) has been coined and fundamentals of a modern bosonization technique 
have been formulated, this model and its various modifications remain at the 
focus of interest in research in strongly correlated systems. Without diminishing 
the role of standard theoretical methods, such as diagrammatic techniques or 
renormalization group approach, one can say that bosonization methods make 
the most powerful theoretical tool for strongly correlated ID systems. All these 
methods, and in particularly the 'canonical' operator bosonization, where the 
creation and annihilation operators of electrons are explicitly represented in 
terms of Bose operators and a 4-fermionic Hamiltonian is eventually diagonal- 
izcd in the bosonic representation, are described by Dmitrii Maslov in lecture 
notes published in this volume. 

The standard operator bosonization is one of the most elegant methods 
developed in theoretical physics. However, by its very formulation it seems 
both limited to and specific for one-dimensional physics. A subject of this 
seminar is to demonstrate the existence and usefulness of an alternative way 
to bosonizc ID interacting electrons, called the 'functional bosonization'. It is 
based on the Hubbard-Stratonovich decoupling of the four-fermion interaction - 
a typical way to "bosonize" a fermionic system in higher-dimensional problems. 
This method was elaborated in different ways in a set of papers [7ll51IUl [TUllll| . 
In this seminar, we will describe such a functional method in the form similar 
to that developed earlier JO., JLL ^^r the treatment of the pure LL as well as a 
single-impurity problem in the Luttinger model. However, we will employ here 
the Keldysh technique (see for reviews ^21) rather than the Matsubara one used 
in 10, 11^. 

The essence of the method is to eliminate a mixed fermion-boson term in 
the action (resulted from the Hubbard-Stratonovich decoupling) by a gauge 
transformation. Such a procedure is exact for the pure ID Luttinger model and 
gives a convenient starting point for including a single backscattering impurity. 

The problem of a single impurity in the LL has been actively investigated 
by many authors J^S, 14, 15, 16, 17, ,18, 19, 20,. One of the main results of these 
considerations [1311141 [THlll9ll20| was the suppression at low temperatures of the 
local density of states (TDoS) at the impurity site and at a distance from the 
impurity JT, and the related suppression of the conductance 13, 14^ and the X 
ray edge singularity ^1 119L OO] . Another prominent result was the dependence 
of the Friedel oscillations ^| El E] on the distance from the impurity. 

In this seminar, we will briefly outline only the most important results for 
the one-impurity problem in the Tomonaga-Luttinger model, while giving a 
slightly more detailed description of the functional bosonization in the Keldysh 
technique that was not previously published. For simplicity, we will only ad- 
dress a single-mode Tomonaga-Luttinger model, with one species of right- and 
left-moving electrons, thus omitting spin indices and considering eventually the 
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simplest linearized model of a single-valley parabolic electron band. We will 
also skip over a physical introduction, referring the reader to lecture notes by 
Dmitrii Maslov published in this volume. 

2 Functional integral representation 

Within the limitations outlined above, the most generic Hamiltonian of inter- 
acting ID electrons in the presence of an external scattering potential v{x) can 
be written as 



where Vq{x — x') is a bare electron-electron interaction. The observable quan- 
tities to be calculated with this Hamiltonian are the tunneling density of states 
(TDoS) and the current as a linear response to an applied field characterised 
by the vector potential A{x,t). We will give the results for the TDoS at the 
end of this presentation, but will mainly described techniques for calculating 
the current j{x,t). It can be written with the help of the Kubo formula in 
terms of the current-current correlation function thermally averaged ((. . .)q) in 
the equilibrium Gibbs ensemble: 



where the current operator is defined in the standard way (in the units h = 1 
here and elsewhere): 



The Kubo formula can be rewritten in Keldysh techniques by defining the con- 
tour Ck which runs from — oo to the observation time t along the upper bank 
of the cut along the real time axis in the complex time plane and then returns 
to — oo along the lower cut: 



Here we have introduced the chronological operator time-ordering the opera- 
tors in a descending order along the contour Ck (with all the times on the lower 
cut considered as later times as compared to those on the upper cut). Although 
the current operators are bilinear in Fermi-operators, for the future usage we 
assume in the standard way that the Fermi-operators anticommute under Tk, 




(1) 




j{x,t) = -—iJ^x,f:)d,;ip{x,t) +h.c. 



(3) 
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assuming additionally that at equal times if^^t) is taken at an infinitesimally 
later moment then ip{t). Such an agreement is consistent with the definition 
of the current (Eq|3I)- In what follows, we will choose a more general contour 
(the Keldysh contour), running from minus to plus infinity above the cut and 
returning below the cut. In the operator language, such an extension of the 
contour corresponds to the insertion of extra evolution operators: having been 
not coupled to the observables (the current operator in our case) , such evolution 
operators above and below the cut simply cancel each other. 

Now, any expression written as the chronological operator average can be 
straightforwardly represented as a functional integral over the fields defined on 
the contour of the time ordering. We introduce the action, 

S[i;,i:] = J dt{i4,M~H[^,ij]} , (5) 
Ck 

with the last term being the the normal ordered Hamiltonian with 'tp^ ^ and 
ip —> ip, where ip and ^ are the anticommuting Fermi fields. Then, the the linear 
response current of Eq. takes the form 

j{x, t)=i j At' j dx' l^j{x, t)j(x', t')) A(x', t') - "^Aix, t), (6) 

Ck 

where the brackets stand for the functional integral 

(...)= J V4,vi^{...)e^^^^^\ (7) 

and the current field is defined by 

— 6 — 

t) = 7, — - ^{x, t)dx^{x, t) + c.c. 

3 Effective action for the Tomonaga - Luttinger 
Model 

For free electrons, the plane waves basis is natural. Assuming that the scattering 
by impurities and electron-electron interaction involve the energy scales much 
smaller than the Fermi energy, the plane wave basis presents a natural starting 
point. Therefore, we may assume that the main contribution to the functional 
integral above comes from the fields representing separately right-moving and 
left-moving electrons 

^{x, t) « t) e'^^" + ^l{x, t) e-*P^" (8) 

where '4>r,l are smooth on the scale. Such a separation is the essence 
of the Tomonaga-Luttinger model, and corresponds to the linearization of the 
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initially parabolic electron band. The right- and left-moving electrons can be 
transformed one to another due to backscattering processes. We shall neglect 
such processes due to the electron-electron interaction, thus keeping only the 
small momentum transfer part of interaction, V{q <C ^pp). This part of inter- 
action is non-trivial by itself as it breaks the Fermi Liquid theory in ID. On the 
other hand we will keep only backscattering in the impurity potential v{x) since 
small-momentum elastic scattering does not result in any qualitative change in 
the Luttinger Liquid behavior. 

Now we make the substitution |(SJ) neglecting higher order derivatives of 
smooth functions and discarding integrals over fast oscillating terms. After 
some straightforward manipulations we come to the action for the Tomonaga- 
Luttinger model: 

Stl — So + Sint- (9) 

The first term describes free electrons in the presence of the external scattering 
potential (which can be also changing in time): 

5„^/d.d.*t(,,,)(^^|^.^^ f^;*))*(.,t). (10) 

We have assumed (here and below) that all the time integrations are performed 
along the Keldysh contour, and introduced the following notations: 

* = (^jfj , *^ = (V^^, , dn/L = dt± vpd, . (11) 

The second term in Eq. Q gives the interaction part of the action 



v'dtn{x,t)Vo{x~ x')n{x' ,t) , (^n = ^-tfj . (12) 

4 The bosonized action for free electrons 

Before dealing with interacting electrons, we convert the action H1U|) for free 
electrons in the presence of the impurity potential v into the action in terms of 
bosonic fields. 

Expanding Eq. (|10|l in powers of v, we obtain the partition function by 
integrating over the fermion fields with the help of the Wick theorem: 

^0 = X! (^^1)2 d^-zcTz' detgL{zi,z'j)detgR{zl,Zj)Y[vizk)v{zl). (13) 

n=0 k=l 

We have introduced the notation z = (x, t) and defined the Green functions of 
left- and right- moving electrons in the standard way (suppressing indices L, R): 

g ^ -i (V(z)V5(z')) = . (14) 
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The matrix structure is defined by the position of the time arguments of ip 
and "0: in and g both of them are, respectively, on the upper or lower 
branches of the contour, and in g^{g^) the first (second) argument is on the 
upper branch while the second (first) is on the lower. We should remind that 
the functional average Q, invoked in the definition H14f) . automatically arranges 
for the time ordering along the Keldysh contour. The Green function obeys the 
equation 



idrjgrjiz, z') =8{x~ x')6{t,t'), 



(15) 



Here S{t,t') = S{t — t') when both arguments are on the same branch of the 
contour and S{t,t') = otherwise; 9,, is defined in Eq. Hll|l . Explicitly 



T 



1 



"^^F sinhvrT 



{t - t'-iO) T 



(16) 



while g^ is different by the sign of iO. The components of g obey the usual 
relation g++ + g = + g^. In what follows, we will perform the standard 
Keldysh rotation, reducing all the appropriate matrices to the triangular form 
with the Keldysh component {g^ = 5^+5^) in the upper right corner, and the 
retarded and advanced components {g'' = (5°)* = g^^ — g^) on the main diag- 
onal (see, e.g., ref. Jl])- We shall assume that the time arguments belonging 
to the upper and lower branch of the contour have, respectively, positive and 
negative infinitesimal shift into the complex plane. Note finally that in the equi- 
librium case presented here the (Fourier transform of the) Keldysh component 
is related to the (Fourier transforms of the) retarded and advanced ones via the 
Fermi distribution function fe{T) as follows: 



g^ie)^{l-2MT)){g^e)^g^ie)) 



(17) 



The same relations are valid for other Green's functions to be considered so 
that, wherever this docs not involve an ambiguity, we will give only the re- 
tarded/advanced components in the explicit form implying that the expression 
like Eq. (|17|l is valid for the Keldysh component. 

To calculate the partition function we use the Cauchy identity |21j : 



det 



1 



sinh(zi — z') 



sinh(zi — Zj) sinh(Zj' — zj) 

(-1)"+^ . (18) 

J_|smh(zi - Zj) 



It reduces Eq. (|13|l for the partition function to the following one: 



.k=l 



Il,<jS{zt-Zj)s{z'i-z'j) 



(19) 
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with 

s{z - z') = sinhTrTU - <' - ^— ^ ) sinhTrTU - t' + \ (20) 

V VF J V VF J 

Introducing the bosonic Green function 

iGo(z-z') = -lns(z-z') (21) 
whose retarded and advanced components are Fourier transform of 

<'^(g,c.) = ^i^^, u;±^Lo±zO, (22) 

one can write 

.2n n ri 



^"^Exni^i^) / nd^^dz^«(zfc)i;(z;) 

ri=0 ^ '-^ ^ ^ / 

xexpp^[Go(z^-z,) + Go(z,'-4)] + ^ ^ Go(z. - z^ )] . (23) 

This is a partition function of the Coulomb gas with the logarithmic interaction 
and it can be represented as the functional integral over the bosonic field ip{x, t): 

Zo = Jvtp exp ^^iSo[ip] +aj dz[ve-'"f + c.c] | . (24) 

Here the free bosonic action So[lp] is defined in terms of the Green's function 
Go of Eq. inj as follows 

5oM = ^ [dzdz'^{z)G^\z - z')^{z') , (25) 



The constant a in Eq. (|24|l absorbs an ill-defined value of Go(x, t) at x = 0,t — 0. 
We use the ultra-violet cutoff which corresponds to the scale of order ef- 



211^ 2ttvf 

Thus we have cast the original free fermion problem into that of interacting 
bosons, represented by the partition function H24|l . The interaction between 
bosons, i.e. the second term in the exponent in Eq. H24() . comes from the back- 
scattering impurity term in the original fermionic problem. The Gaussian action 
for noninteracting bosons, Eq. (|25|l . can be explicitly written via x and t in the 
standard form 



HirvF 



{dt^f -v%{d,,^f , (26) 



although Eq. (|25|l is no less convenient, especially for the generalization for the 
interaction. The main advantage of the bosonization, either in the standard or 
functional form, is that including the quadric electron-electron interaction does 
not substantially change the free action. 
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5 Gauging out the interaction 



The first step in including the interaction term H12() into the bosonization scheme 
is to perform the Hubbard-Stratonovich transformation which can be symboh- 
cally written as 

exp{i5int[V^,V']} = J I?0exp|-^0VV + '/'*^*|- (27) 

Note that the auxihary 'Hubbard-Stratonovich' (HS) bosonic field </> here is 
different from the field (p in Eqs. (|24l) and H26|l . Substituting this representation 
into the full action (jS)), we bring the partition function to the following form: 

Z = J X>(/)P*tx)^e*'^o[0]+iS[*,0] ^ (•28) 
The action S'[\E', for fermions interacting with the HS field is given by 

si*,*i./d.*.,,,('^«-;* *(-')• M 

To cast this integral into the form identical to that of the previous section we 
apply the local gauge transformation, 

ipr,{z) = ipr,{x, t) ip^ix, t) e'^(^^*) with idr,0r,{x, t) = (t>{x, t) , (30) 

which removes the bosonic field from the diagonal part of the action H29|l 
but at a cost: the off-diagonal terms are rotated with the factors e**^, and the 
Jacobian of the transformation J changes the quadratic in 4> part of the action. 

It is shown in Appendix A that the Jacobian J of the gauge transformation 
l|5n|l can be represented as 

In J[cf)] = -\j dzdz>(z) H(z, z') (f>{z') . (31) 

The polarization operator H is given in the random phase approximation (RPA) 

by 

H= ^ H^, n^iz - z') = ig^{z - z')g^{z' - z) , (32) 

ri=R,L 

where gr,{x — x' ,t — t') is the free electron Green function given by Eq. IjlGfl. It 
is well known that the RPA is exact for the Luttinger Liquid Jj. Note that we 
give in Appendix A a very simple and straightforward proof of this. 

The quadratic in cj) contribution of the Jacobian, Eq. H31|l . should be added 
to the quadratic term in Eq. (|27|l . This results in the free bosonic action with 
the kernel corresponding to the screened interaction: 

S[(f>]=-^ l'dzdz'(j){z)V-\z- z')(j){z'), V-^^Vq-^+H. (33) 
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The above expression should be understood in the operator sense: V and 11 are 
the operators whose kernels are defined with the appropriate time and spatial 
dependence on the Keldysh contour. 

Using the expressions for the retarded and advanced components of g^j (with 
e± = £ ± iO), 

dR^il' ^) = (£± ~ ""Fqy^ , //"(q, = (e± + VFq)~^ , (34) 

one finds the appropriate components of the polarisation operator, 

nr/a, . 1 Vpq „r/a, ^ 1 Vpq 



2TT uj± ^ Vpq 2tt uj± + Vpq 

and thus the total polarisation operator as 

W^%q,.) = -l 7^^^ ^ , (35) 

Assuming that the Fourier transform of the forward-scattering pair interaction 
only weakly depends on momentum, i.e. Vo{q^2p^) « const = V, and substi- 
tuting Eq. H35(l into the free bosonic action H33|l . one finds the components of 
the free HS bosonic propagator as follows: 

,,2 „2 



(36) 

V^(<7,^) = tanh (^) [V%,Lo) - V^iq,L.)] . 

Here we introduced the renormalized velocity v which defines the effective cou- 
pling constant g: 

V =vp^ , 5 = — . (37) 

TT V 

Therefore, the gauge transformation H3U|I reduces the action in H28|) to 

S = S[9]+S[¥,^-,e]. (38) 
Its fermionic part is given by 

^[*t, ^.e^^ j d,$t(,) ^ ^_9r^ « ^ ^^^^ ^3g^ 

with 9 = 9p( — 9^, while its bosonic part S[9] is defined via the field (j) by Eqs. 
H33|) , H3()|l and (|37|l . It is convenient to write it explicitly as an integral over the 
field 9, which is straightforward since 9 is linearly related to the field 4> as in 
Eq. H30() . Thus we arrive at the following explicit expression for S[9] in Eq. (|38() : 

S[9] = ^J dzdz'9{z)Gg\z - z')9{z') . (40) 
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The Gaussian kernel of this interaction, Gb, can be represented as 



Gb — G — Go 



(41) 



where Go is defined by Eqs. (|21|) and l|22|) . while G has the standard trian- 
gular matrix structure in the Keldysh space, with the Fourier transform of its 
retarded/advanced component given by 



(42) 



i.e. it differs from Gq only by substituting v for vp in the denominator. Finally, 
the expressions for fermion density and current in new variables become 



2,71 



27r 



(43) 



which can be seen by keeping source terms coupled to original fcrmionic fields 
when calculating the Jacobian. 

The effective action (I38II is quadratic in fcrmionic fields which can now be 
easily integrated out. Before doing so, let us stress that the representation 
of Eqs. H38|l - H4()|l seems to be more convenient for some problems than the 
fully bosonized action. To perform the fcrmionic integration, we note that the 
fcrmionic part of the action, Eq. H39|l . differs from that for the free electrons, 
Eq. H10() . only by the substitution v — > ve~^^. Therefore, repeating the same 
procedure as in the previous section, we represent this part of the action with the 
help of the bosonic field ip so that the full action in Eq. H38|l goes (in symbolical 
notations) to 



s[ip,e]^^eGs'e + ^^G^'ip + a 



Introducing the new bosonic field, Q = 9 + ip, and noting again the relation 14111 
we arrive at the standard fully bosonized action, 

S'[e] = i fdzdz'e{z)G-\z- z')e{z') + aldz\^viz)e-'^^''^ +c.c^ , (44) 



with Eq. H43I) transforming to 



n = no + -;^dxQ, j ^ --^dtQ . 



(45) 



6 Tunnelling density of states near a single im- 
purity 

As an application of the formalism developed above we will calculate the density 
of states in the vicinity of a single impurity in the Luttinger Liquid characterised 
by the following local time-independent potential: 

vix) = Xvf6{x) (46) 
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where A is the dimensionless impurity strength. The potential (|46ll should be 
substituted into the action (|44|l . Then one integrates out the fields with x ^ Q 
which results in the local action in terms of Q{t) = Q{x = Q,t): 

Simp = ^ y dMt'e(t)G,r;^p(t - t')e{t') + 2iaXvF J dt cose(i) (47) 

Ck Ck 

where the Fourier transform of the retarded/advanced components of the inverse 
Gaussian kernel, G'imp(i), are obtained by integrating the Green's functions of 
Eq. (|42|l over all momenta: 

^r/a, X f 4:nvF 2tt 

We now employ the self-consistent harmonic approximation (see, e.g., |22l 
\2'6\). i.e. substitute the impurity cosO term with the quadratic one: 

i2avF\ J dt cose(t) -^A J dte^(t) 

The coefficient A is to be found from the condition that this substitution is 
optimal. 



d_ 
OA 



2avFX{cose) - ^A{e^) 



0, (49) 



where the averages are taken with the effective action symbolically represented 
as 



'S'cff 



ie(Gr4 + A)e. (50) 



Solving self-consistently Eq. (|49|l with the action H5Uf) (which involves preserving 
the proper analytical structure, with A being of the standard matrix structure 
in the Keldysh space) one finds with the logarithmic accuracy: 

Ar/^^±eFl^^X^"^ ^±eFX^. (51) 

Including the source terms of Eq. H45|l does not cause any principal difficulties 
but makes the transformations more cumbersome. Therefore, we do not describe 
such transformations in the framework of this seminar presentation. Instead, 
we simply present the results for the tunnelling density of states [following from 
the long time asymptotics of the full electron Green function calculated in the 
SCHA with the bosonized action H38|) ] as a function of the distance from the 
impurity. We refer the reader interested in detail to our previous publication 
[TT] (albeit for the Matsubara rather than the Keldysh functional integrals). 

We present the explicit expressions for the TDoS smoothed over the length 
scale much larger than in three different regions: 
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A 



-1 



-1 



X 



£9 



A^i < J < £"^ 



min(a;, A 



(a) 
(&) 



(52) 



The regions of different behavior of the tunnehng density of states are sketched 
in the figure. There and in Eq. 152|l . i = e = s/ef, and the renormal- 

ized impurity strength A = lA*"'"! is given by Eq. 1)51(1 . Equation ((52b ) describes 
TDoS in the vicinity of impurity, in full correspondence with the original re- 
sults of Kane and Fisher .13 obtained for the TDoS at a; = 0, i.e. exactly 
at the impurity. In addition, we have established here the TDoS dependence 
on the impurity strength A = A^^. The region of applicability of Eq. ((52b .') 
corresponds to the diagonally hatched region in Fig. 1. Equation ((521 1 gives 
the TDoS at very large distances from the impurity. As expected, it coincides 
with a well-known result for the TDoS in the homogenous Luttinger liquid. Its 
region of applicability is horizontally hatched in Fig. 1. In the intermediate 
region, vertically hatched in Fig. 1, the TDoS depends both on the energy and 
the distance from the impurity. This analytic dependence given by Eq. I(52b) 
describes the crossover from the impurity-induced dip in the TDoS to the bulk 
behavior. Finally, the unhatched region for e < a corresponds to small energies, 
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e < T, where the energy dependence saturates (by e ^ T) in aU the three hnes 
of Eq. (PI). 

In conclusion, we have demonstrated how to develop the formalism of bosoniza- 
tion based on the functional integral representation of observable quantities 
within the Keldysh formalism. We have derived in this way the fully bosonized 
action for the interacting electrons in the presence of the scattering potential 
given by Eq. (|44|l . and illustrated its usage on the example of the TDoS on 
a single impurity, Eq. (|52|) . Let us stress finally that the intermediate repre- 
sentation of Eqs. 13811 14U|) . which still contains the part quadratic in fermion 
fields appears to be more convenient for some problems than the fully bosonized 
action. 

Acknowledgements. We gratefully acknowledge support by the EPSRC 
grant GR/R95432. 



A Jacobian of the gauge transformation 



The Jacobian of the gauge transformation H30() can be defined as J = JrJl 
with 



Trln[l-ig„^ 



(53) 



where the Green functions of non-interacting right- or left-moving electrons, 
obeying 9^(7^ = /, are given by Eq. IjlGI) . Note that in the matrix components 
g^{t) and g^{t) the time argument should be understood, respectively, as tzLiO. 
The exponent in Eq. H53I) can be represented as infinite series in the HS-field 0: 



2 

In J,, [(/}]= -"^—T^t: [9r,<f>r 



(54) 



The n-th order term in is proportional to the loop with n external lines 
corresponding to (^'s, each loop being built of the n Green functions gr/'- 



/- n 
Y[dz,(j){z^) 
-i=l 

where the n-th order vertex is given by 

n 

Tn {zi, Zn) = JJ^ gri {Zi', Zi+i) , 



Tn (^1 7 ■ ■ • 1 -^n) 7 



(55) 



(z„+i = Zl) 



Introducing the new variables ^ = e'^'^'^'"'^ *J , one represents the Green func- 



tions of Eq. (|16|l as 



^2 
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so that the vertex becomes 



r„ (6; oc 7„ (^i; ...;^„)TT^j, 7«==TT7 — 7 — , (56) 

1=1 1=1 ^ ' 

with the boundary condition Cn+i = Ci- Since only the symmetric part of the 
vertex contributes to the integral H55|l . we may symmetrize 7^: 

where is an absolutely antisymmetric polynomial of the N-th order which 
depends on n variables. A simple power counting gives N = n{n — 3) /2 while 
the possible minimal order of such a polynomial is Nmin — n{n + l)/2. This is 
not self-contradictory for n = 1 and n = 2 only. All other terms must therefore 
be equal to zero. The term with n = 1 is cancelled due to electroneutrality. The 
only non-vanishing vertex then is the one with two legs: 

In Jn M = ^ y"dzdz' 0(z)n^(z; z')^(z'), (58) 

where 

J\r,{z] z') ^ i9n{z - z')g^{z' - z) . (59) 
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